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ABSTRACT 

We analyze and interpret SOHO/MDI data on oscillation frequency changes 
between 1996 and 2004 focusing on differences between activity minimum and 
maximum of solar cycle 23. We study only the behavior of the centroid frequen¬ 
cies, which reflect changes averaged over spherical surfaces. Both the f-mode and 
p-mode frequencies are correlated with general measures of the sun’s magnetic 
activity. However, the physics behind each of the two correlations is quite differ¬ 
ent. We show that for the f-modes the dominant cause of the frequency increase 
is the dynamical effect of the rising magnetic field. The relevant rise must occur 
in subphotospheric layers reaching to some 0.5 - 0.7 kG at a depth of about 5 Mm. 
However, the implied constraints also require the field change in the atmosphere 
to be so small that it has only a tiny dynamical effect on p-mode frequencies. For 
p-modes, the most plausible explanation of the frequency increase is a less than 
2% decrease in the radial component of the turbulent velocity in the outer layers. 

Lower velocity implies a lower efficiency of the convective transport, hence lower 
temperature, which also contributes to the p-mode frequency increase. 

Subject headings: Sun : Helioseismology, solar variability 

1. Introduction all of solar cycle 23. Information about 

cycle-dependent changes in solar oscilla¬ 
tions includes data on the mean multiplet 


We now have data on the evolution of 
solar oscillation frequencies covering nearly 
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frequencies, z/£„, the multiplet structures 
described by the ak,en coefficients, corre¬ 
sponding mode amplitudes and widths. In 
this paper, we focus on changes in the v. 

The correlation between p-mode fre¬ 
quencies and the magnetic activity cycle 
was hrst observed during the declining 
phase of the cycle 21 (Woodard & Noyes, 
1985) and conhrmed during the next cycle 
by a number of independent studies (Lib- 
brecht & Woodard 1990; Woodard et ah 
1991; Bachmann & Brown 1993; Elsworth 
et al. 1994; Regulo et ah 1994; Chap¬ 
lin et al. 1998). The increase of f-mode 
frequencies with rising activity was discov¬ 
ered during the rising phase of the current 
cycle. 

The physical origin of oscillation fre¬ 
quency increases with rising activity has 
been a matter of controversy. The first ex¬ 
planation, given by Goldreich et al. (1991, 
hereafter GMWK), was that the dominant 
cause of frequency growth with activity 
is the effect of an averaged small scale 
magnetic held changing the frequencies di¬ 
rectly through the perturbed Lorentz force 
and indirectly through the induced pres¬ 
sure change. An objection to this expla¬ 
nation was raised by Kuhn (2000) who ar¬ 
gued that direct measurements of the rms 
held in the sun’s photosphere (Lin, 1995; 
Lin & Rimmele, 1999) preclude the held 
growth required in this picture. Instead, 
he proposed that the main ehect causing 
p-mode frequency rise is a decrease in tur¬ 
bulent velocity due to the rising held’s in¬ 
hibition of convection, which is believed to 
be the main ehect of the magnetic held on 
the sun’s interior structure (e.g. Spruit, 
2000 ). 

As for the f-modes, Antia et al. (2000) 


hrst noted that the frequency rises with ris¬ 
ing activity, and that it is roughly propor¬ 
tional to the frequency itself. Such a be¬ 
havior could be accounted for by a solar 
radius decrease. The number quoted by 
these authors was 5 km during the three 
years of the rising phase of the cycle. This 
result was broadly conhrmed by Dziem- 
bowski, Goode & Schou (2001, hereafter 
DGS) who analyzed SOHO/MDI data and 
found that the f-mode frequency shifts may 
be explained by two components: one be¬ 
ing similar to that found by Antia et al. 
(2000), and the other growing more rapidly 
with frequency. DGS pointed out that 
the former component cannot arise from 
a shrinking of the photospheric radius, but 
rather from a layer located between 4 and 
8 Mm beneath the photosphere. They also 
suggested that the shrinking is caused by 
an increase of the radial component of a 
small-scale magnetic held beneath 8 Mm 
below the surface having a 1 -10 kG level 
rms value. 

The result presented by DGS was criti¬ 
cized by Antia (2003) who argued that the 
f-mode signal found in SOHO/MDI could 
be completely accounted for by an annual 
variation of a non-solar origin. We do not 
agree with his criticism. However, we still 
regarded it useful to reconsider our inter¬ 
pretation having now a much larger data 
set in hand. Indeed, the interpretation pre¬ 
sented in this paper is diherent: we now at¬ 
tribute the observed rise of f-mode frequen¬ 
cies directly to the rise of the magnetic held 
in the layers sampled by f -modes - that is, 
the outer 8 Mm. 
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2. Frequency changes between 1996 
and 2004 from SOHO/MDI data 

Libbrecht & Woodard (1990), who first 
determined activity related p-mode fre¬ 
quency shift for modes over a broad range 
of degrees, £, noted that most of the fre¬ 
quency dependence of the shift is described 
by the inverse of the mode inertia, 
which they called mode mass. We thus ex¬ 
press the frequency shifts in the form 

Az/fe = (1) 

where lin is dimensionless mode inertia, 
which is calculated assuming a common 
normalization of the radial mean displace¬ 
ment in the photosphere. Such calcula¬ 
tions require a solar model. In this work, 
we use model S of Christensen-Dalsgaard 
et al.(1996). The adopted normalization 
is such that /io,i 9 = 1- Values of I^n de¬ 
crease with i. The n-dependence is more 
complicated. At low n, there is a sharp in¬ 
crease. The minimum inertia is reached at 
intermediate orders, which at £ = 10 corre¬ 
sponds to n = 22 and a frequency of about 
3.8 mHz. The p-mode data extend up to 
i = 200 and cover a frequency, u, range of 
1.1 — 4.5 mHz. For these latter modes, the 
mode dependence is essentially reduced to 
a simple //-dependence (see Fig. 4 in DSG). 
The lack of a separate ^-dependence tells us 
that the sources of frequency shift must be 
localized in the outer layers above the lower 
turning point for the modes in the sample, 
or at least for the modes that matter. 

We emphasize that we treat the f-modes 
separately, because even in the outer lay¬ 
ers these modes have vastly different prop¬ 
erties than those of p-modes at the same 
frequency. Hence, we cannot expect the 


same 7 (z/) dependence for both types of 
modes. The kernels for calculating q’s re¬ 
sulting from changes in the magnetic held, 
turbulent pressure, and temperature cal¬ 
culated by Dziembowski & Goode (2004, 
hereafter DG) are indeed very different 
for these two types of modes. In partic¬ 
ular, for p-modes the dominant terms in 
the kernels are proportional to |V ■ 
where ^ denotes the displacement eigenvec¬ 
tor, while for f-modes the dominant term 
is ^l^l/r ^ IV ■ ^|. In the next section, we 
will, as in DGS, consider representations 
including the “radius” for f-modes, but we 
will argue against a signihcant role for it. 
Both types of 7 (z/)-dependence are helio- 
seismic probes of the averaged changes over 
spherical surfaces in the subphotospheric 
layers during the activity cycle. However, 
they are independent probes. 

The plots in Fig.l show the frequency 
averaged 7 ’s for all available datasets 
from SOHO/MDI measurements calcu¬ 
lated from frequency differences relative 
to the hrst set from activity minimum of 
cycle 23. For comparison, in the bottom 
panel we show the monthly sunspot num¬ 
ber taken from the National Geophysical 
Data Genter . For the f-modes, we £t con¬ 
stant (//-independent) values for 7 . For 
the p-modes, we fit a three term Legendre 
polynomial series. We used all available 
frequencies in the two averaged sets. The 
frequency differences are weighted by the 
inverse of the sum of the squared errors. 
The similarity in the behavior of the p- 
and f-modes seen in the two upper pan¬ 
els might suggest that the source of the 
changes is the same in both cases but, as 
we shall see, this is not true. 
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3. The 7 (z/)-dependence for f- and p- 
modes 

The z/-dependence yields a clue to the 
physics of frequency change. By averaging 
frequencies over hve sets covering nearly 
one year of the solar minimum phase and 
ten sets covering two years of the maxi¬ 
mum phase, we averaged out the annual 
changes that are apparently non-solar in 
origin (DGS; Antia 2003). Most of the 
short time scale frequency changes are a re¬ 
flection of variations in solar activity, but 
our aim here is to explain the dominant 
source of the change between minimum 
and maximum. The price for averaging 
the p-mode data is a that is about a 
factor of three larger than those for indi¬ 
vidual sets. Note in Fig.l that the dis¬ 
persion among the averaged sets is much 
larger than the errors. 

Figure 2 shows individual 7 £„ values 
with the la error bars and the Legendre 
polynomial hts. The £t depends on trun¬ 
cation order, N^j., of the polynomial, but 
the results stabilize at Atr = 3 for the f- 
modes, and Atr = 7 for p-modes. The ro¬ 
bust feature of the 7 (z/) dependence for the 
f-modes is the gradual decrease of 7 be¬ 
tween u = 1.37 and 1.74 mHz, correspond¬ 
ing to the Grange of 185 — 300. 

Including an Ai?/-term representing 
variations of the f-mode radius does not 
lead to stable results. The radius Rf ~ 
O. 99 i? 0 , defined in DGS, corresponds to 
the bottom of the layer where all the f- 
modes in the data set are trapped. The 
numbers in Table 1 demonstrate the erratic 
behavior of Ai?/ and 7 with increasing Atr. 
Within the error bars, the result obtained 
with Atr = 0 is the same as quoted by 


DGS. With more terms, the fit is improved 
but the inferred values become meaning¬ 
less. Low and stable 7 values are ob¬ 
tained only after excluding the Ai?j term. 
We thus abandon the idea that the rise of 
f-mode frequencies is caused by a shrink¬ 
ing of the radius beneath the bottom of 
the f-mode zone. With the data available 
in 2000 , when the first interpretations of 
the f-mode frequency changes were given, 
the data did not allow for more than a two 
parameter fit. Furthermore, we found that 
averaging the frequency data taken during 
solar minimum and maximum sets was es¬ 
sential for the present inference. In Section 
4, we will argue that the dominant part of 
the frequency increase is due to an in situ 
rise of the magnetic held. 

Table 1: Fitting the f-mode frequency 
shifts to 


Au 


3ARf 

2 ^~ 


iy + 


/ 


using the Legendre polynomial series trun¬ 
cated at Atr for xi^)- 


Atr 

Ai?/[km 

7 


0 

0 

0.70± 0.02 

3.18 

1 

0 

1.19± 0.04 

1.56 

2 

0 

0.91± 0.06 

1.30 

3 

0 

1.06± 0.08 

1.25 

0 

-4.79 

0.41± 0.04 

2.14 

1 

5.41 

1.94± 0.15 

1.39 

2 

-7.90 

-0.50± 0.58 

1.26 

3 

10.23 

3.06± 2.60 

1.26 


The robust feature of the 7 (z/)-dependence 
for the p-modes is the steady increase be¬ 
yond u = 2 mHz. We stress that the sig- 
nihcant decreasing trend in 7 (z/) over the 
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1.4 to 1.74 mHz frequency range found 
for the f-mode cannot be replicated if one 
looks for a common 7 (z/) dependence for 
all the modes, which naturally would be 
dominated by the more abundant p-modes. 
For both p- and f-modes, higher frequency 
means a stronger sampling of the outer¬ 
most layers. Therefore, the opposing be¬ 
havior of the two types of modes at the 
high frequency end of the spectrum sug¬ 
gests that different physical effects are re¬ 
sponsible for the frequency increase corre¬ 
lated with rising solar activity. 

The smallest values of for the p- 
modes are signihcantly higher than those 
for the f-modes. Undonbtedly, the rel¬ 
atively higher is caused by temporal 
fluctuations in activity dnring the maxi- 
mnm phase (for the f-modes snch fluctu¬ 
ations are closer to being within the rel¬ 
atively larger error bars). However, the 
relatively high for the maximum phase 
with respect to the minimum might also be 
due to an inadeqnacy of the £t in which 
the mode dependence in 7 ^^ comes only 
throngh freqnency. For instance, an ad¬ 
ditional mode dependence is expected dne 
to changes bnried in the deep layers lo¬ 
cated below or near the turning of certain 
p-modes in the sample. To check, we plot 
the residuals against the position of the 
mode turning points, as determined by 

„ , „ , ImHz 

= (£ + 0.5)-, 

V 

and show the resnlts in Fig. 3. We see 
that for fly > 50, corresponding to a turn¬ 
ing point of depth of less than 40 Mm, 
the residuals are on average somewhat less 
than zero, while they are greater than zero 
for fly < 50. This means that a de¬ 
tectable contribntion to the p-mode fre¬ 


qnency changes arises in the layers reach¬ 
ing down to a depth of 40Mm. However, 
most of the contribution arises in much 
shallower layers. Note that there is no vis¬ 
ible contribntion from the vicinity of the 
bottom of the convective envelope. 

There are various potential contribntors 
to the frequency changes described by the 
7 (z/) fnnctions. These include the mechan¬ 
ical effect of the spatially averaged changes 
of the magnetic held, as well as the effects 
of such changes on the convective and ther¬ 
mal strnctnre of the outer layers. These 
were discussed, e.g., by GMWK and DGS. 
A full description of all these contribntors 
requires five unknown functions of depth. 
We note that for the parts of the interior 
most robustly probed, the 7 (z/) from ob¬ 
servations are not snfhciently accurate to 
even think about a formal inversion. All 
that may be done is to fit simple fnnctional 
dependences for each specified contributor. 
This is what we will do here following ap¬ 
proaches adopted by GMWK and DGS in 
their assessments of the held required to 
explain p-mode frequency changes. 

4. Variational expressions for 7 

Hamilton’s variational principle is our 
tool for interpreting the freqnency changes. 
It was employed, for instance, by GMWK 
and by DG, whose integral formulae link¬ 
ing the 7 ’s to changes in the magnetic 
held, turbulent pressure and temperature 
are used in the present paper. The ap¬ 
proach adopts an adiabatic approximation 
for oscillations, which may not be fully jus¬ 
tified in the part of the onter layers of 
our interest. DG also provided expressions 
for calculating frequency splittings repre¬ 
sented by the even-a coefficients. Here we 
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give a summary of only the formulae that 
are relevant for calculating the centroid fre¬ 
quencies. 

The underlying variational expression 
for the angular frequency, u}{= 27rz/), shift 
used here is 

A, _ A (Dp Dm + Dy) 

2Ioj ■ ^ ^ 

The quantities Dp, Dm-, and D^ repre¬ 
sent contributions of pressure, magnetic 
field, and turbulent velocity, respectively, 
to mode frequency. The quantity 

1 = j £xp\^\^ = R^pi, (3) 

is the mode inertia and ^ is the Lagrangian 
displacement vector calculated in a stan¬ 
dard spherical model. Here we have ig¬ 
nored terms resulting from changes in 
gravity and large scale velocity fields, since 
both were found to be negligible. 

Only the rms values averaged over 
spherical surfaces contribute to the changes 
in the mean frequencies, 9. The pressure 
term may be written in the form 

Dp = J d3a;p[S+(r-l)|div||2], (4) 

where S = and denotes the co¬ 

variant derivative. The Lagrangian change 
of Dp, which is calculated with S{pd^x) = 
0, requires evaluation of the pressure and 
density perturbations, Sp and Sp respec¬ 
tively, induced by changes in the magnetic 
helds and correlated changes in the tur¬ 
bulent velocity. Such changes arise from 
the magnetic field’s inhibiting of convec¬ 
tion. For the spherically symmetrical part 
of perturbation considered here, the den¬ 
sity perturbation is not determined by the 


condition of mechanical equilibrium. Here, 
following DG, we express Sp in terms of 6p 
and the Lagrangian temperature perturba¬ 
tion, ST. We treat changes in the mag¬ 
netic field, turbulent velocity, and temper¬ 
ature as independent sources of frequency 
changes, even though they are physically 
linked. Modeling the effect of magnetic 
held on convection is still not well under¬ 
stood, and our approach is to use simple 
reliable physics to derive helioseismic con¬ 
straints on advanced models. 

Both the magnetic and velocity helds 
are treated as being statistically random 
with the net ehect on radial structure re¬ 
sulting from their square averaged com¬ 
ponents. The vertical component was al¬ 
lowed to be diherent from the two hori¬ 
zontal components. Thus, the covariance 
matrix for the magnetic held is written in 
the following form, 

= Sij[SjrM^{r) + 

-M^{r){Sjg + Sj^)], (5) 

where Sij is the Kronecker symbol. We will 
use an unsubscripted S to denote the La¬ 
grangian changes in solar depth dependent 
parameters, while A refers to changes in 
the global parameters. An analogous form 
was adopted for the turbulent velocities. 

pv~v] = p5ij[5jrTi^ (r) + 

+ ( 6 ) 

Upon assuming mechanical equilibrium in 
a thin layer, DG (see also GMWK) calcu¬ 
lated the change in pressure, Sp, resulting 
from changes in the random magnetic held 
{SM^, SAi^) and in the random velocities 
{ST^, SV^). The evaluation of the den¬ 
sity change requires consideration of the 
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thermal balance. DG calculated ADp as¬ 
suming an isothermal response of the layer, 
and separately evaluated the contribution 
to the 7 ’s from the temperature change, 
6T. 

The total dynamical effect of the mag¬ 
netic held change on frequencies consists of 
ADp calculated with ST = 0 and 


ADm 


-A 


d^x 


An 

2div|*B ■ {B ■ 

+ Idiv^P 




-\B\ 


(7) 


These two terms combined in Eq. 2 lead 
to the following expression for 7 , 


7m 



d. 


phot 


1 Mm 


K 


V 

M,0 


fSM^\ 

\l kGV 


-H 


+^M,0 


fsM^y 

\1 kGV. 


/r Hz, ( 8 ) 


where dphot denotes the depth beneath the 
photosphere. The explicit expressions for 
and in terms of radial eigen¬ 

functions of the mode are given in Eq.(67) 
of DG. For the case of f-modes we have 
77 m 0 ~ t^M,o > 0 - For p-modes we have 
77 m 0 ^ |77moI- Fig. 4 shows examples of 
the kernels. 


^M,0 


— 2 (77m,0 + 277m, oJ 




for calculating 7 due to isotropic changes 
in the held. 

Gomparing the kernels for the two f- 
modes shown in Fig. 4 with the behavior 
of the 7 (zz) shown in the upper panel of 
Fig. 2, we conclude that, if the rise of the 
magnetic held is responsible for the f-mode 


7 ’s, which decrease between u = 1.28 and 
1.74 mHz, the growth must occur predom¬ 
inantly beneath 2.5 Mm, where the kernel 
of the lower frequency mode {i = 160) has 
higher value. 

The kernels for the p-modes plotted in 
the lower panel of Fig.4 were calculated for 
7 = 10. However, in these outermost layers 
the kernels of all p-modes in the sample are 
virtually ^independent. The two kernels 
have very low values at depths beneath 2.5 
Mm. It is thus clear that the observed p- 
and f-mode frequency increases cannot be 
simultaneously explained by an increase in 
the magnetic held. This conclusion does 
not depend on the assumed isotropy of the 
held. Isotropy would have been an essen¬ 
tial assumption if the held were the cause 
of the p-mode frequency rise. 

The plots in Fig. 5 illustrate the large 
diherences between p- and f-modes kernels 
at the same frequency. We see that the f- 
modes in the MDI sample are far more sen¬ 
sitive to magnetic held changes in the outer 
few Mm beneath the photosphere than the 
p-modes. We should also note a signih- 
cant diherence between pi and the higher 
order p-modes above dphot = 1 Mm. This 
diherence, however, has virtually no con¬ 
sequence because data on pi and even p 2 
modes are irrelevant for probing this out¬ 
ermost layer. 

The overall dynamical ehect of the tur¬ 
bulent velocity changes is calculated in a 
similar way to that for the magnetic held. 
To the induced change of Dp calculated as¬ 
suming ST = 0, we add the change in the 
velocity term. 


AD^ = -A 


d^xp\{v-V)^\^ 


When these two terms are used in 


• (9) 

Eq. 2, 
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we get, after integration over spherical sur¬ 
faces, 


Iv 


d 


d 


phot 


+K. 


H 


1 Mm 


K, 


V 


1 km^s ^ 


km^s 

jj, Hz. (10) 


Again, we do not give expressions for 
and here. They were given in Eq.(59) 
of DG. In that paper, we plotted the ker¬ 
nels and with k = 0 referring to 
the centroid changes and A; > 1 to the even- 
a coefficients. Unfortunately, there was a 
numerical mistake in those plots. Here we 
also corrected the error in the previous cal¬ 
culations, reducing the values by a factor 
of about 9. In Figure 6 , we plot corrected 
kernels for 7 due to changes in the vertical 
component of the velocity for the same four 
modes that were selected for use in Fig. 4. 
The absolute values of the are smaller 
(factor i for the f-modes, and much less 
than that for the p-modes). Comparing 
the plots the for f- and p-modes, we con¬ 
clude that only in the latter case may one 
expect a signihcant effect from changes in 
the turbulent velocities. The effect should 
arise mainly above the depth of 1 mM, 
where the highest velocities are expected. 
Large effects of turbulence on solar f-modes 
found in a number of investigations (e.g. 
Murawski, 2000) concerned modes of much 
higher degrees than considered in this pa¬ 
per. 

A decrease in the vertical component of 
the turbulent velocity remains the most vi¬ 
able explanation of the dominant part of 
the p-mode frequency increase correlated 
with the magnetic activity cycle, because 
high frequency modes preferentially sam¬ 
ple the layers where we expect the largest 


changes in the turbulent velocity. 

The decrease in the velocities also 
means a decrease in the efficiency of con¬ 
vective energy transport, hence a decrease 
of temperature in the outer layers. 

As GMWK hrst observed, an isobaric 
increase of temperature causes an increase 
in p-mode frequencies. The same is true 
for f-modes, but there the effect is much 
smaller. We express the 7 ’s due to the La- 
grangian change in temperature, 6T, in the 
form 


again referring readers to DG for an ex¬ 
plicit expression for Kt- The plots of Kt 
for the same four modes as in Fig. 6 are 
given in Fig. 7. Here, we also decreased 
the values by a factor of about 9 relative 
to corresponding plots in DG. The ker¬ 
nels’ behavior is similar to that seen in 
Fig. 6 . In the shallow subphotospheric lay¬ 
ers, where we may expect the largest rela¬ 
tive changes of temperature caused by de¬ 
creased efficiency of convective transport, 
only p-modes have substantial amplitudes, 
and it is only for these modes that a tem¬ 
perature decrease must be considered to be 
a potentially important contributor to the 
frequency increase. 

5. Changes in snbphotospheric mag¬ 
netic field from changes in the 
f-mode frequencies 

The dynamical effect of the rise of the 
magnetic held seems to be the only pos¬ 
sible explanation for the observed f-mode 
frequency increases. Formally, one may ex¬ 
plain the behavior of the 7 (zz)’s seen in 



Fig. 2 in terms of a decrease in the turbu¬ 
lent velocity, but the decrease would have 
to be nearly constant with depth, and the 
amount would have to be unacceptably 
large at depths greater than say, 2 Mm. 

Probing the depth dependence of the 
magnetic held based on f-modes has mod¬ 
est precision. Modes in the [180, 300] l- 
range, for which we have a signihcant de¬ 
termination of 7 , effectively sample the 
magnetic held down to a depth of only a 
few (5-6) Mm, but the probing precision is 
not high because of the spread in the indi¬ 
vidual frequency shifts. As a hrst attempt, 
we seek the isotropic held change, = 

in the form of a truncated 
power series of the depth below the tem¬ 
perature minimum, dmin = dphot + 0.476 
Mm, 


at dint = 4 Mm and j > 20 (M2 in Fig.8). 
However, similar values of were reached 
at higher dint and lower j. One such exam¬ 
ple (M3) is shown in Fig. 8. 

We see that a concave shape of 7 (z/) is 
reproduced only with the model allowing 
d(i?^) < 0, but we do not regard this hnd- 
ing to be signihcant. Rather, we would 
blame some inadequacy in our model of 
the small-scale held in the atmosphere. 
What we regard to be signihcant is that 
the f-mode frequency increase between so¬ 
lar minimum and maximum requires both 
an average held increase of some 0.5 - 0.7 
kG at a depth of about 5 Mm, and a much 
smaller increase close to the photosphere. 

6. The dominant source p-mode 
variations 


( 12 ) 

k=0 

Including terms up to A; = 2 enabled 
us to reach similar values of to those 
of the three-term Legendre polynomial ht 
of the 7 (z/) dependence. The resulting 
d(R^)(dmin) dependence does not look re¬ 
alistic because we found 5(5^) < 0 near 
the photosphere (see curves denoted Ml in 
Fig. 8), as though the near-surface held 
decreases with rising activity. With this in 
mind, we tried a form of the 6(drain) 
function forcing d{B'^) > 0 everywhere. 
We chose 


5{B^) 


6(B )int if drain ^ dint 

l5(S^)lnt (iSf)'’ it dmln < <ilnt 


(13) 

with adjustable parameters d(i?^)int, d^t, 
and j. The lowest y^ of 1.69 was reached 


Changes in the magnetic helds inferred 
from f-mode data have only a very small 
ehect on p-mode frequencies. This is illus¬ 
trated in Fig. 9 where in the upper panel 
we show the 7th-order Legendre polyno¬ 
mial hts of the measured frequency diher- 
ence after removing the ehect of the held 
according to model M2 (here the choice of 
the model is not important for our conclu¬ 
sions). Only in the lower frequency regime, 
below V = 2 mHz, are the increases in 
the 7 ’s with decreasing u reduced. This 
suggests that the averaged dynamical ef¬ 
fect of the magnetic held rises at a depth 
of a few Mm, and is responsible for an 
appreciable part of the frequency increase 
of low frequency p-modes. However, we 
should stress that, as we may see in Fig. 
2, the signihcance of 7 (z/) in this part of 
the p-mode spectrum is questionable. In 
any case, most of the p-mode frequency in¬ 
crease with rising activity requires a diher- 
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ent explanation. 

The high frequency part of the 7 (z/) de¬ 
pendence, which is really signihcant, may 
be explained only by invoking an effect act¬ 
ing preferentially close to the sun’s photo¬ 
sphere. The dynamical effect of the grow¬ 
ing magnetic held is excluded by measure¬ 
ments of the averaged photospheric held 
and by the f-mode data. What remains 
to be considered is an inhibiting ehect of 
the held on convection leading to a lower 
turbulent velocity and temperature in the 
outermost layers. The two ehects are ex¬ 
pected to be signihcant only close to the 
photosphere. The question to answer is 
how much of a reduction is required to ac¬ 
count for the observed frequency changes. 

We hrst consider the ehect of lowering 
the turbulent velocity. In fact, only the 
vertical component of it matters because 
the horizontal components hardly ahect p- 
mode frequencies. The squared averaged 
vertical velocity in Eq. (10) is represented 
in the form of a truncated power series of 
the depth beneath the temperature mini¬ 
mum, dmin, 

= (14) 

k=0 

It turned out that it suffices to include 
terms up to /c = 2 to ht the measured y’s 
with a better than that from the sev¬ 
enth order polynomial 7 (zz). The hts are 
compared in the upper panel of Fig. 9. 

In the lower panel of Fig. 9, we show 
the inferred d(T'^)(dphot) dependence with 
the 1(7 error bars from the least square £t. 
The changes required to account for the p- 
mode frequency rise are naturally higher 
than those assessed by DG, as the ker¬ 
nels they used were grossly exaggerated 


(by about a factor of nine due to an error, 
which has been fixed here), but the crude 
estimate made therein results in an error 
that is smaller than should have been ex¬ 
pected. With our present corrected and 
precise analysis, we get higher numbers, 
but they are not unreasonably high. The 
comparison with the model values shows 
that what is required is less than a 2.5 
% decrease in T^. That is less than a 
1.3% decrease in the rms vertical compo¬ 
nent of the turbulent velocity. We do not 
believe this number is in conflict with ob¬ 
servations. 

A potentially more difficult problem 
arises from implication of the reduced con¬ 
vective efficiency on the effective tempera¬ 
ture. According to a crude estimate, based 
on mixing-length approximation and an 
Eddington atmosphere given by DG, a 1 
percent decrease in the convective velocity 
is associated with a relative temperature 
decrease ranging from 1 x 10“^ at dphot = 1 
to 3 X 10“^ at dphot = 0 Mm. The values 
are about one-half of what is needed to 
account for the p-mode frequency increase 
by the temperature effect alone. Thus, the 
required velocity reduction is smaller. As¬ 
suming a 0.65% percent reduction in the 
r.m.s. turbulent velocity and adopting the 
relation (^D(O), we find 8 K for 

the required decrease in the effective tem¬ 
perature, which seems unacceptably high. 
Nonetheless, we believe that the inhibiting 
effect of the magnetic held on convection 
is the cause of the p-mode frequency in¬ 
crease correlated with increasing activity. 
We cannot conceive of a more plausible ex¬ 
planation, and we blame the problem re¬ 
garding the effective temperature change 
on inadequacies in our treatment of energy 
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transport in the convective zone and in the 
atmosphere. 

7. Conclusions 

We analyzed all available SOHO/MDI 
data to study the behavior of the mean so¬ 
lar frequencies with varying magnetic ac¬ 
tivity. Averaged over their respective fre¬ 
quency ranges, the time variations of p- 
and f-mode frequencies show the same pat¬ 
tern. Both are strictly correlated with 
sunspot number. The difference is seen 
when the mode dependence of the fre¬ 
quency shifts between activity maximum 
and minimum is compared. The quantities 
we compare are the shifts scaled by mode 
inertia, that is the 7 ’s. There is a slight 
residual mode dependence for p-modes in¬ 
dicating that there is a contribution to the 
shift arising in deeper layers, but still well 
above the bottom of the convective enve¬ 
lope. In the frequency ranges where 7 (zz) 
is well determined, the two types of modes 
exhibit opposite trends with increasing fre¬ 
quency: growing 7 ’s for p-modes and de¬ 
clining ones for f-modes. We determined 
different scenarios as the explanation of the 
dominant source of the frequency changes 
in these two cases. 

We considered two possible sources of 
the mean frequency changes: ( 1 ) dynami¬ 
cal effects of the changing, average, small- 
scale magnetic held; ( 2 ) effects of turbulent 
velocity and subphotospheric temperature 
changes caused by the impeding effect of 
the held on convection. In our analysis, 
we relied on formalism developed by DG. 
We also corrected a numerical error in es¬ 
timates presented in that work. 

We demonstrated that the main part of 
the f-mode frequency shifts is explained by 


the growth of the subphotospheric mag¬ 
netic held. The relevant growth takes place 
in the layers at depths of 2.5-5 Mm rising 
to about 0.5—0.7 kG. The detailed implica¬ 
tions regarding both shallower and deeper 
layers are uncertain beyond the sharp de¬ 
crease in the held required toward the sur¬ 
face. Formally, the best ht is obtained if 
there is a slight decrease in the mean pho- 
tospheric held in the outermost layer with 
increasing activity, but we do not believe 
that this is realistic beyond saying the held 
growth there is small. 

Because of its location, the held caus¬ 
ing f-mode frequency rise has only a minor 
ehect on p-modes. The weak held rise in 
the outermost layers is also consistent with 
the direct measurements of the mean pho- 
tospheric held (Lin 1995, and Lin and Rim- 
mele 1999). This outermost layer is where 
the dominant source of the frequency p- 
mode change resides. Thus, we can exclude 
the held as the source of p-mode changes. 
Attributing the p-mode frequency shifts to 
a decrease in turbulence, we found that 
this requires less than a two percent de¬ 
crease in the rms value of the vertical com¬ 
ponent, as calculated in a model of solar 
convection. Impeded convective hows also 
cause lower temperatures in the outermost 
convective layers. Lowering the tempera¬ 
ture causes a frequency rise because the 
ehect of cooling is more than compensated 
by the resulting contraction. 

The difficulty of the proposed explana¬ 
tion of the p-mode frequency shifts comes 
from the implied decrease of the ehective 
temperature by some 8 K between solar 
minimum and maximum phase. Since our 
estimate is based on a very crude model, 
we do not regard this problem as essential. 
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but we think it calls for a closer study with 
advanced models of the convective zone 
and atmosphere of the sun. 

Let us return to the question posed 
in the title of the DGS paper; Does the 
sun shrink with increasing magnetic activ¬ 
ity ? The answer following from our present 
analysis is: yes, it does but not because 
of the changes at depths beneath 8 Mm, 
which were previously suggested to lead to 
shrinking with the rate of about 1.5 km/y 
during the rising activity phase. The data 
covering the whole high activity phase al¬ 
low a multiple parameter £t of the fre¬ 
quency change dependence on frequency, 
Az/(z/). We found in Section 3 that there 
is no stable solution for the shrinking rate 
at the depth of 8Mm. The dominant effect 
influencing mass distribution in the outer¬ 
most layers is the decrease of the turbulent 
pressure and temperature with increasing 
activity, and both effects cause shrinking. 
The implied shrinking amounts to about 1 
km between solar minimum and maximum. 
The roughly 1% decrease of the squared 
turbulent velocity in the outer 1 Mm below 
the photosphere must be compensated by a 
0.1 percent density rise because turbulence 
contributes about 10 percent of the total 
pressure. We would obtain the same esti¬ 
mate of the shrinking by attributing part of 
the p-mode frequency increase to the tem¬ 
perature decrease. 

This estimate for the radius change is 
anti-correlated with activity, and has the 
opposite sign and is much smaller than 
the change of the photospheric radius de¬ 
rived by Emilio et al. (2000) from the di¬ 
rect measurements based on SOHO/MDI 
intensity data. They determined 5.9±0.7 
km/y for the rate of the photospheric ra¬ 


dius increase during the rise phase of cy¬ 
cle 23. We stress , as in DGS, that the 
two rates are not directly comparable be¬ 
cause the change we derive refers to radius 
at constant mass and not to constant op¬ 
tical depth. 

A cooler and smaller active sun, whose 
increased irradiance is totally due to ac¬ 
tivity induced corrugation, has been advo¬ 
cated for years by Spruit (e.g. 1991, 2000). 
Our results support his picture. 
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Fig. 1.— The values of the averaged 
7 ’s, which are global helioseismic mea¬ 
sure of solar activity, are derived from 
38 SOHO/MDI data sets compared with 
monthly sunspot numbers shown in the 
bottom panel. Note that the p-mode 7 ’s 
(mid panel) closely replicate changes in 
the sunspot nnmber (lower panel) and also 
other general measnres of solar activity 
during the cycle. The behavior of the f- 
modes (npper panel) is similar, but the val¬ 
ues are less significant. The larger errors 
are mainly a consequence of the fact that 
the f-mode spectrum is an order of magni¬ 
tude sparser. 


This 2-column preprint was prepared with the 
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Fig. 2.— Frequency dependence of the 
7 ’s derived from the frequency difference 
between averaged frequencies from solar 
maximum phase (2000.4 - 2002.4) and the 
minimum phase (1996.3 - 1997.3). The 
lines represent fits using truncated Legen¬ 
dre polynomial series. Subscripts at y’s 
denote the order at which the series was 
truncated. The quoted values of s-re 
calculated per degree of freedom. 
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Fig. 4.— Kernels for calculating the y’s 
arising from isotropic changes in the mag¬ 
netic field according to Eq. 8, for two se¬ 
lected f-modes (upper panel) and p-modes 
(lower panel), and at two selected frequen¬ 
cies, are plotted as fnnctions of depth be¬ 
neath the photosphere in the outer part of 
the standard solar models. 
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Fig. 3.— Residuals of for p-modes after 
removing the 77 (z/) fit (see Fig. 2) plotted 
against position of the lower turning point. 
The value (£-l-0.5)/z/ = 50 corresponds to a 
depth of 40 Mm. The vertical line indicates 
the bottom of the convective envelope. 



Fig. 5.— Similar to Fig. 4 but for modes 
of nearly same frequency of u = 1.74 mHz. 
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Fig. 6 .— Kernels for calculating 7 ’s aris¬ 
ing from the change in the random velocity 
field according to eq. 10 for two selected 
f-modes (upper panel) and p-modes (lower 
panel) and at two selected frequencies plot¬ 
ted as functions of depth in the outer part 
of the standard solar models. 
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Fig. 7.— Kernels for calculating the fre¬ 
quency shifts due to temperature increase 
for p-modes at selected frequencies (lower 
panel), and f-mode modes at selected de¬ 
grees (upper panel). 
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Fig. 8 .— The upper panel shows the 
depth dependence of the mean square held 
change between solar maximum and max¬ 
imum for models Ml, M2, and M3. The 
errors in individual values of for the 

Ml model are large, growing from 0.03 at 
dphot = 0 to 0.3 (kG)^ at dphot = 5.5 Mm. 
In the lower panel, the bolder lines (in red) 
show the f-mode q’s calculated for each of 
the three models. The thinner lines (in 
blue) show 7 ( 1 /) functions obtained by ht- 
ting truncated Legendre polynomial series 
as described in Section 3. 
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Fig. 9.— Solid and dashed black lines in 
the upper panel show 7 (z/) functions ob¬ 
tained by htting the Legendre polynomial 
dependence to measured p-mode frequency 
shifts before and after removing the ef¬ 
fect of magnetic held changes according to 
model M2 (see upper panel of Fig. 8 ). 
Solid red triangles denote individual values 
of 'Jin calculated assuming that the p-mode 
frequency shift is caused by a decrease in 
the turbulent velocity. The lower panel 
shows the inferred values of the change in 
the mean squared turbulent velocity and 
compares them with the 1 percent decrease 
of those values calculated in a model of 
the solar convective zone of Abbett et al. 
(1997). 
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